(1)
摘 要：生态系统模型模拟是研究全球变化背景下区域气候变化对生态系统影响的重要方法之一，大多数气候模型预测全球大部分区域的NPP将增加，然而未来气候情景下近百年的长时间尺度下NPP及其稳定性的变化却鲜有研究。本数据结合建立在生理生态学基础上的CEVSA模型，通过强迫法发展了遥感驱动的生态系统过程模型——CEVSA-RS模型，采用的基于区域气候模型第4版（RegCM4.6）和CMIP5的HadGEM2-ES数据情景，输入驱动的气候变量，包括时、空间分辨率分别为0.25°，3 h的气温、降水量、云量和空气相对湿度，后将其处理为以旬为时间步长，0.1°空间分辨率，使用Chinacover2010的土地覆盖数据，模型输出为空间分辨率为0.1°的2006—2099年的NPP及其变化趋势，以及全期（2006—2099年）、前期（2006—2035年）、中期（2036—2065年）和远期（2066—2099年）NPP稳定性多年均值。通过共享本数据集，以期为积极开展减缓和适应气候变化行动，如双碳行动，提供具有科学和现实意义的数据参考。
关键词：净初级生产力NPP；CEVSA-RS模型；未来气候情景；稳定性
DOI: https://doi.org/10.3974/geodp.2023.02.05
CSTR: https://cstr.escience.org.cn/CSTR:20146.14.2023.02.05
数据可用性声明：本文关联实体数据集已在《全球变化数据仓储电子杂志（中英文）》出版，可获取：
https://doi.org/10.3974/geodb.2023.06.01.V1 或 https://cstr.escience.org.cn/CSTR:20146.11.2023.06.01.V1.
植被净初级生产力（Vegetation Net Primary Productivity，NPP）表征大气碳的固定和生物质的累积[1]，是研究全球陆地生态系统对气候变化响应的基础[2]。驱动碳汇变化的气候因素具有较强的时空异质性[3]，生产力的波动会对碳固定也有一定的影响[4]。如何应对全球变化下，升高的气温，变化的降水格局，判断NPP的稳定程度，是目前的热门研究课题之一。
生态系统模型模拟是研究区域生态系统对气候变化的响应的重要方法之一[5, 6]，通常分为气候统计模型、遥感参数模型和生态过程模型[7, 8]。以小尺度精细试验分析为基础的生态过程模型，从机理层面出发，具有较好的未来变化预测能力，通过结合遥感参数模型，实现尺度融合，以追求更高精度的生态系统碳动态模拟[9—12]。将水文、大气、二氧化碳、土地覆盖、物种组成和氮沉降等过程考虑进气候碳循环模型中，可以提高模型模拟近期和未来气候变化的能力，是未来模型的发展方向之一[13]。遥感驱动的生态系统过程模型（CEVSA-RS），对生态系统生产力的模拟具有较高的精度[14, 15]。本数据通过使用关闭遥感驱动的模型，在气候驱动下模拟未来气候情景陆地生态系统碳循环时空动态，为碳循环模拟研究提供了方法基础。
模型输入选择CMIP5中基于不同温室气体排放情景，其中，RCP4.5模拟经济均衡发展模式，作为气候变化的“底线”；RCP8.5作为在不实施有效气候变化减缓措施条件下，代表本世纪最恶劣的气候情景。
本数据仅选择了这两种情景数据，作为模型输入，应用CEVSA-RS模型模拟了未来气候情景下2006—2099年全国植被NPP，分析了NPP的大小、变化趋势和稳定性，以期为全国陆地生态系统的碳汇管理，生态恢复措施实施及生态系统管理策略制定等提供科学基础。目前鲜有对NPP稳定性的定义，本数据通过变异系数的概念，将NPP稳定性定义为变异系数的倒数，由此表征NPP面对气象因子等变化因素的稳定程度。同理，目前也少有NPP稳定性的数据，以此数据为后人提供参考，为缓解气候变化的政策提供基础。
《未来百年气候情景下中国植被净初级生产力稳定性估算数据集》[16]的名称、作者、地理区域、数据年代、时间分辨率、空间分辨率、数据集组成、数据出版与共享服务平台、数据共享政策等信息见表1。
| 条目 | 描述 |
|---|---|
| 数据集名称 | 未来百年气候情景下中国植被净初级生产力稳定性估算数据集 |
| 数据集短名 | RCPsNPPChina |
| 作者信息 | 陈惺，中国科学院地理科学与资源研究所 中国科学院大学，chenx.20s@igsnrr.ac.cn 王军邦，中国科学院地理科学与资源研究所 中国科学院大学，jbwang@igsnrr.ac.cn 何启凡，九江学院旅游与地理学院，2051936579@qq.com 王春雨，齐鲁师范学院，1871302580@qq.com 叶辉，九江学院旅游与地理学院，fever2cn.huiye@outlook.com |
| 地理区域 | 3°51′N—53°33′N，73°33′E—135°05′E |
| 数据年代 | 2006—2099年 |
| 空间分辨率 | 0.1° |
| 数据格式 | .tif |
| 数据量 | 4.26 MB（压缩后） |
| 数据集组成 | 一共16个栅格数据文件，包括RCP4.5和RCP8.5情景下多年平均的NPP及其变化趋势、NPP稳定性的多年平均栅格数据 |
| 基金项目 | 国家自然科学基金（31861143015, 31971507）；中国科学院—青海省人民政府三江源国家公园联合研究专项（LHZX-2020-07） |
| 出版与共享服务平台 | 全球变化科学研究数据出版系统 http://www.geodoi.ac.cn |
| 地址 | 北京市朝阳区大屯路甲11号100101，中国科学院地理科学与资源研究所 |
| 数据共享政策 | 全球变化科学研究数据出版系统的“数据”包括元数据（中英文）、通过《全球变化数据仓储电子杂志（中英文）》发表的实体数据集和通过《全球变化数据学报（中英文）》发表的数据论文。其共享政策如下：（1）“数据”以最便利的方式通过互联网系统免费向全社会开放，用户免费浏览、免费下载；（2）最终用户使用“数据”需要按照引用格式在参考文献或适当的位置标注数据来源；（3）增值服务用户或以任何形式散发和传播（包括通过计算机服务器）“数据”的用户需要与《全球变化数据学报（中英文）》编辑部签署书面协议，获得许可；（4）摘取“数据”中的部分记录创作新数据的作者需要遵循10%引用原则，即从本数据集中摘取的数据记录少于新数据集总记录量的10%，同时需要对摘取的数据记录标注数据来源[17] |
| 数据和论文检索系统 | DOI，CSTR，Crossref，DCI，CSCD，CNKI，SciEngine，WDS/ISC，GEOSS |
本数据的模型输入使用了2006—2099年中国陆地生态系统气候情景数据[18, 19]。该数据是基于区域气候模型第4版（RegCM 4.6）和CMIP5中的HadGEM2-ES数据情景，本数据使用其中的中（RCP4.5）、高（RCP8.5）排放气候情景数据。原数据的时、空间分辨率分别为0.25°，3 h，包括气温、降水量、云量和空气相对湿度，后将其处理为以旬为时间步长，0.1°空间分辨率，用于模型输入。
本数据使用的土地利用类型输入数据是基于Landsat TM/ETM、HJ-1卫星数据、野外调查数据得到的2010年中国土地覆盖遥感数据，原数据空间分辨率为30 m，全国尺度一级分类精度为94%，二级分类精度为86%[20, 21]。后对其二级分类进行归并，并采用邻近值法空间重采样得到空间分辨率为0.1°的模型输入数据。详细的模型输入输出见表2。
本工作基于中、高排放气候情景下2006—2099年全国陆地生态系统气候情景数据，采用2010年中国土地覆盖遥感数据，根据发展的基于卫星遥感驱动的CEVSA-RS模型模拟计算了不同未来气候情景下的NPP数据。
| 数据类型 | 指标 | 时间分辨率 | 空间分辨率 | 单位 |
|---|---|---|---|---|
| 输入气象数据 | 气温（Tas） | 旬均值 | 0.1° | ℃ |
| 降水（Prc） | 旬均值 | 0.1° | mm | |
| 相对湿度（Hum） | 旬均值 | 0.1° | % | |
| 云量（Clo） | 旬均值 | 0.1° | % | |
| 输入土地利用数据 | 土地覆盖数据（Chinacover2010） | 年值（2010年） | 0.1° | 无 |
| 输出变量 | 植被净初级生产力（NPP） | 年值 | 0.1° | g C·m?2·a?1 |
| 输出变量 | NPP稳定性 | 年值 | 0.1° | 无量纲 |
将得到的NPP数据进行趋势分析和变异型分析，得到对应的全时段NPP变化趋势和分时段（前期、中期、远期）的NPP稳定性空间数据，数据处理具体包括以下几个部分：
（1）趋势分析能更好地表达NPP对气候变化的响应，公式如下：
(1)
式中，Slope为像元尺度线性回归方程的斜率，Xi为第i年的NPP值；n为研究的时间长度。
（2）稳定性可以表示一个系统保持不变或者规律性变动的稳定程度，变异系数本身表达波动性[20]，因此通常采用变异系数倒数的绝对值量化稳定性，计算公式如下：
(2)
式中，SX、MX和STDX分别表示给定时间范围内生态系统生产力、气候系统中气温和降水等量的稳定性、均值和标准偏差，其中，气温数据有负值，因此通过取绝对值的方法以消除气温的负值影响。i代表第i年，j代表第j个像元点。
本数据集一共16个栅格数据文件，包括RCP4.5和RCP8.5情景下多年平均的NPP及其变化趋势、NPP稳定性的多年平均栅格数据[22, 23]。以30年为一个时期，分为近期、中期、远期。
（1）RCP45代表RCP4.5气候情景，RCP85代表RCP8.5气候情景；
（2）npp_mean代表2006—2099年的NPP均值；
（3）npptrend20062099代表2006—2099年像元级别的NPP变化趋势；
（4）P20062099代表2006—2099年像元级的NPP变化趋势的显著性水平；
（5）R220062099代表2006—2099年像元级的NPP变化趋势的决定系数；
（6）stable代表稳定性多年均值，其中20062035代表2006—2035年近期的多年平均NPP稳定性，20362065代表2036—2065年中期的多年平均NPP稳定性，20652099代表2066—2099年远期的多年平均NPP稳定性，20062099代表2006—2099年全时段的多年平均NPP稳定性。
关于单位的说明：
（1）NPP：g C·m?2·a?1；
（2）NPP稳定性：无量纲；
（3）NPP变化趋势：g C·m?2·a?1；
（4）显著性水平、决定系数：无量纲。
本数据在像素尺度准确量化了两个未来气候情景下的NPP的空间异质性及其变化趋势以及NPP稳定性的多年均值。2006—2099年中国陆地生态系统植被净初级生产力多年均值总体上呈现东高西低，南高北低的空间格局（图1）。其中，在RCP4.5的情景下，全国的NPP总值达到了4.41 Pg C·a?1，相比于RCP8.5情景，RCP4.5情景略高，NPP的空间差异较大，相同气候区内的NPP较相似，如青藏高原区的NPP总量为0.57 Pg C·a?1，是中国四大气候区中总NPP量最少的区域，热带-亚热带季风区是四大区域中NPP总量最大的，为2.09 Pg C·a?1。
图1 RCP4.5情景下2006—2099年中国陆地生态系统NPP均值及其变化趋势图
下载（注：该图基于自然资源部标准地图服务网站审图号为GS(2016)1594号的标准地图制作，底图无修改）
然而，根据NPP的变化趋势空间格局，除了青藏高原气候区，其余气候区大部分地区都为下降趋势（图1）。未来百年的NPP变化不容乐观。全国NPP总量的94年年际变化趋势结果显示，RCP4.5情景和RCP8.5情景下都呈现显著下降的趋势（图2），其中，中等排放情景的下降速率要大于高排放情景。并且，值得注意的是，在21世纪60年代前后，NPP总量的变化将出现转折，转折点的形成原因值得探讨。
图2 2006-2099年中国陆地生态系统NPP总量年际变化图
下载分时段的NPP稳定性空间分布格局结果显示（图3），2006—2099年NPP稳定性青藏高原区呈现高值，而温带季风区，温带大陆气候区呈现低值；2006—2035年（近期）的NPP稳定性呈现南高北低；2036—2065年（中期）NPP稳定性同样呈现南高北低，但南部的高值相对近期较高；2066—2099年（远期）青藏高原区NPP稳定性相对高，出现大量的NPP稳定性较差的区域。其中，值得注意的是，青藏高原区呈现出较强的空间异质性。
图3 RCP4.5情景下不同时段的中国陆地生态系统NPP稳定性均值图
下载该模型在发展过程中基于中国陆地生态系统通量观测研究网络（ChinaFLUX）观测数据进行了验证评价研究。模拟结果与内蒙古、长白山和禹城站点的涡度相关碳通量观测数据对比表明，模型能够解释总初级生产力（GPP）观测值季节变化的58%—94%[14]；同时与MODIS的GPP遥感产品（MYD17A2H）表现出了较好的一致性[14]，说明CEVSA-RS模型具有较高的可靠性[16]。
本数据不仅包括CEVSA-RS模型模拟的两种未来气候情景下NPP空间格局，还包括近期、中期、远期NPP稳定性的空间格局。可用于评估不同分区下的NPP总量，还可用于评估NPP应对气候变化的稳定程度。CEVSA-RS模型已应用于多区域，多种生产力的评估，数据结果可靠。本数据研究气候变化情景下植被生产力及其稳定性的变化，可为考虑区域差异开展诸如生态保护与修复及双碳行动等减缓和适应气候变化的举措等提供科学参考。根据数据结果显示，未来百年的NPP总体趋势表现为显著下降，并且中远期的NPP稳定性空间分布结果显示，我们应该重视未来百年的NPP变化，从对NPP总量的定量分析转化为NPP稳定性的定量分析，并且青藏高原地区展示了独特的区域异质性，我们应该有差异地决策不同气候区对未来气候变化如何进行响应。
作者分工：王军邦对数据集的开发做了总体设计；陈惺采集和处理了NPP稳定性数据；叶辉和何启凡采集和处理了NPP数据；王军邦设计了模型和算法；王春雨做了数据验证；陈惺，Watson, A. E.撰写了数据论文。
利益冲突声明：本研究不存在研究者以及与公开研究成果有关的利益冲突。